Kondo effect in an antiferromagnetic metal 
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We study the fate of a spin- 1/2 impurity in the itinerant antiferromagnetic metallic phase via 
a renormalization group analysis and a variational calculation. The local moment - conduction 
electron interaction hamiltonian in an antiferromagnetic metal is spin non-conserving. We show 
that for a general location of the impurity, the Kondo singularities still occur, but the ground state 
has a partially unscreened moment. We calculate the magnitude of this residual moment and the 
variation of the spin polarization with energy for a substitutional impurity as a function of the 
staggered magnetization. The usual Kondo effect only occurs if the impurity is placed at points 
where the magnetization is zero. 

PACS numbers: 



I. INTRODUCTION. 

The dual nature of electron behavior in antiferro- 
magnetic (AFM) heavy fermion (HF) materials is not 
well understood. On one hand, in many intermetallic 
compounds containing elements with partially filled f- 
orbitals, the heavy mass originates from the collective 
screening of the local moments by the conduction band 
via the Kondo effect (for a review see RefsQIii). Emer- 
gence of magnetic order, on the other hand, requires in- 
teraction between the unscreened moments. 

In many materials exhibiting coexistence of antiferro- 
magnetism and heavy mass the RKKY interaction be- 
tween the local moments is comparable to the Kondo 
temperature4ii^, and the AFM order involves only part of 
the local moments at low temperatures^. This suggests 
a ground state where a fraction of the full unscreened 
moment appears as ordered staggered moment, while the 
rest continues to be compensated by the conduction elec- 
trons via the Kondo effect^. In materials with several 
/-electrons per unit cell, such as several U-containing 
compounds, the coexistence of AFM and HF behavior 
may be understood as originating from Kondo screen- 
ing of one /-species and ordering of another. In systems 
with a single /-electron, for example many Ce-containing 
intermetallics, the coexistence of magnetic ordering and 
screening originates with the same local moment. 

Motivated by this picture we consider in this paper 
the Kondo effect for a single impurity coupled to a band 
of conduction electrons that order antiferromagnetically. 
We believe that this is the first step towards develop- 
ing a low-energy theory for the coexistence of the AFM 
and HF behavior. Just as the essential aspects of sin- 
glet formation and mass enhancement in paramagnetic 
heavy fermions follow from the analysis of the screen- 
ing of a local moment in an electron bath, the salient 
features of competition between antiferromagnetism and 
Kondo screening can be understood from the analysis of 
our model. We write down a simple model hamiltonian 
for this problem and show that in a mean-field theory 
for the lattice, our model has precisely the local symme- 
tries required to consider the competition between the 



two phenomena. We find that only partial screening of 
the impurity spin takes place due to spin non-conserving 
interaction vertices, with the unscreened moment frac- 
tion that depends on the amplitude of the AFM order. 

A major advantage of our model is that it allows a 
variety of analytic approaches that make the underlying 
physical picture more transparent. Previous analyses of 
the competition between antiferromagnetism and screen- 
ing for the Kondo lattice or equivalent models have been 
largely numerical. Several authors considered a half-filled 
Kondo lattice with anisotropic Kondo couplingii^. In 
that case the AFM phase is insulating, rather than metal- 
lic, in contrast to the experimental situation. Moreover, 
the coexistence regime does not exist for the isotropic 
Kondo coupling that we consider, although it may ap- 
pear under applied magnetic fields^. A metallic phase 
with coexisting Kondo screening and AFM order was re- 
cently found in variational Monte Carlo calculations for 
the Kondo Lattice Modeli^, however, Kondo screening 
was argued to remain the same in the paramagnetic and 
antiferromagnetic phases, in contrast to our findings be- 
low. 

Our study also has fundamental importance beyond 
the connection to the heavy fermion physics. Our main 
conclusions, (i) that the Kondo screening of an impurity 
placed in an itinerant antiferromagnet (a) depends on 
the location of the impurity within the unit cell; (b) for 
substitutional impurity on a simple lattice is incomplete 
with the magnitude of the residual moment related to the 
amplitude of the AFM order; (ii) that the spin-dependent 
local density of states refiects this incomplete screening 
can all be tested experimentally. It is important to em- 
phasize that, for a system without nesting, after the onset 
of the AFM order and doubling of the unit cell there re- 
mains a large Fermi surface with non-vanishing density 
of states, and one may naively expect a complete Kondo 
screening. The existence of residual unscreened moment 
is due to the spin non-conserving vertices in the magneti- 
cally ordered state that alter the interaction between the 
conduction electron bath and the local moment. 

The remainder of the paper is organized as follows. In 
the next section we introduce our model and discuss its 
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salient features. The structure of the interaction vertices 
and the renormahzation group analysis are presented in 
Sec lIIII We then present the variational ground state 
ansatz in Sec lIVi and follow it by the discussion of the 
density of states on the impurity site. 



II. MODEL. 

Our model consists of an electronic band with an itin- 
erant AFM modulation coupled to an isolated impurity, 
with the hamiltonian H = Hc + Hj + Hk- The conduc- 
tion band without magnetic order is described by 



(1) 



where e (k) is the energy band dispersion, and cj^ ^ is the 
creation operator for an electron eigenstate in a param- 
agnet. The AFM order is imposed via 



H, 



■ym ( 4jCk+Q,T - ^jCk+qa ) + h.c. 



(2) 



where m is the mean field staggered magnetization, cho- 
sen along the spin z-axis, ym is the AFM gap, and Q is 
the ordering wave- vector. 

We consider an isolated impurity with spin S, at po- 
sition Ri coupled antiferromagnetically to the fcrmion 
spins. The Kondo interaction is 



(3) 



where J is the (Kondo) exchange coupling between the 
impurity spin and the conduction electrons, cr is the vec- 
tor of Pauli matrices, and V'K-R) = J2k a exp(ik • R), 
is the creation operator in real space. 

The justification for this model is as follows. Suppose 
we consider the full Kondo lattice Hamiltonian including 
the RKKY interaction and allow for the possibility of 
antiferromagnetic order so that the local spin at each 
site may be written as 

S, = (Sq) exp(iQ • R,) + [S, - (S)q cxp(iQ • R,)] . (4) 

Now consider a mean field theory (for example the dy- 
namical mean- field theory). At the first iteration of the 
self-consistent solution, one must solve the local "impu- 
rity" problem, without making this decomposition at the 
impurity site but making it everywhere else to generate a 
self-consistent bath. Then Sq.® follows as the periodic 
spin-dependent potential on the conduction electrons due 
to the staggered magnetization, 7m oc J{Sq). Therefore 
the local problem is precisely the problem with the sym- 
metries of Eqs.([T|)-(I3]). For a fully self-consistent solution 
we would have to relate 7m to the sub-lattice moment 
generated at the impurity site(s) due to it; this second 
part is not attempted here. The self-consistency affects 
the quantitative details, but not qualitative physics of our 



conclusions. Our finding that there is only partial Kondo 
screening in the local model confirms the consistency of 
the approach. 

In the absence of the Kondo coupling, the Hamilto- 
nian, Eqs.(IT])-([2]) can be readily diagonafized using a Bo- 
goliubov transformation. This introduces four species of 
fermions, described by the band indices n = ±, in addi- 
tion to the two spin projections on the z-axis a = ±1. 



(5) 



with the Bogoliubov factors 



(6) 



Under this transformation the band hamiltonian becomes 
Hb^ Hc + Hi ^ En{k)al^^^anka , (7) 



k,n,Q 



with the energy dispersion E± (k) = ^k ± V'^k + 7^""i^, 
and the shorthand notation ^k = 5 (^(k) + £(k + Q)) 
and (5k = ^ (^(k) — e(k -f- Q)). The momentum summa- 
tion in Eq. ([7]) is over the first magnetic Brillouin Zone 
(MBZ). In the absence of nesting, ^k 7^ 0, the electrons 
remain ungapped, and the MBZ contains a Fermi sur- 
face (FS). Superficially, since the density of states at the 
Fermi level is finite, this may suggest that Kondo screen- 
ing should occur as in a normal metal. We show below, 
however, that this is not the case since the Kondo inter- 
action vertices have nontrivial spin structure. 

The effective hamiltonian is now written in terms of 
the operators a„,k.Q as H = Hh + Hk where Ht is given 
by Eq. ([7]), and the Kondo interaction is 



H 



K 



_J_ 



^ S • CTq/j /„ka/;^,k'/3"lkaan'k'/3 X^^) 



f+ka = Uk - aUkC*^^' 

f-ka = avk + tike"'^' 



(8b) 
(8c) 



Again, the momentum summation in Eq. (|8ap is over the 
first MBZ. Eq. ((8a|) shows that the effective Kondo inter- 
action, written in terms of the cigenstatcs of the antifcr- 
romagnct is, in general, anisotropic, k-dcpendent, as well 
as dependent on the location of the impurity in the lat- 
tice because of the structure factors, fnka ■ The impurity 
spin is coupled to fermions in both bands, and this inter- 
action is weighted by the structure factors. Despite this 
apparent complication, we show below that the invariant 
vertices have a factorizablc form so that the problem is 
still tractable. 
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FIG. 1: a) Particle-Particle channel, (b) Particle-hole chan- 
nel contributions to the vertex and (c) Contribution to the 
Pseudo-Fermion Self Energy 



III. RENORMALIZATION GROUP ANALYSIS 



A. Spin structure of invariant vertices 



To identify the invariant vertices we employ the 
Abrikosov pseudo-fermion representation for the local 
moment, S = J2fj,u 'S'/^f V'/i^f i where S^,^ is the spin ma- 
The integral equation for the interaction ver- 
tex is represented in Fig. [1] and contains contributions 
from the particle-particle (Fig.[TJa)) and the particle- hole 
(Fig. [TJ6)) channels. For Fig. [Ija) we have (fcs = h = 1) 



+ ^ (V' "^w) /"ka/>q„"G„.„. (*c^.,q) V (-^c.« + ^u;)^;^"/y (q,k'), 

a n /3 q cjs 

where G and F are the band fermion and the pseudofermion propagators respectively. The contribution of particle-hole 
channel is obtained by changing the sign of the frequency in one of the fermion propagators. 

We find the solution for the vertex in the factorized form F^f^,^, (k, k') = Taa';Pi3' fnkaf*'^^' ■ The spin matrix 
^aa' ■,f3p' is separated into distinct symmetry channels which couple the impurity moment to the fermions. From 
Eq. © we find 

rL';/3/3' = ^0 {^aa' S^p' ) + ^1 (c^aa' ' ^130' ) + {^2 + ^I) f^L' ^fifl' ~ (^2^ ^ ^s) <„' ^^p' (10) 

I 



where the label T refers to the sum of the particle-hole 
and particle-particle channels. In the equation above, the 
vertices Fi and F5 conserve the total spin, S + cr, of the 
impurity-conduction electron system, while the vertices 
F2, F3, and F4 do not. This spin non-conservation is the 
main reason for the nontrivial results found below. 



B. Impurity location. 

The structure of the vertex depends on the location 
of the impurity. Consider, for example, a system with 
a square lattice, and the AFM ordering wave vector 
Q = (tt, tt). If the impurity is located at the point where 
the magnetization vanishes, midway between two near- 
est neighbors (R, = (1/2, 0) or (0, 1/2)), exp (Q • r/) = i. 
Then f^^af'i^na = 1 and no change in Kondo screening 
occurs on going from the paramagnetic to the AFM state. 
This special behavior takes place only if the impurity is 
located at the point where the magnetization vanishes. 



For a general position of impurity /^ka/nka 7^ so that 
finite spin non-conserving vertices exist. In the following, 
we consider a substitutional impurity, R = 0. 



C. One-loop Renormalization Group Analysis 

We analyze the equations for the vertices at the 
one loop level, replacing the full propagators in Eq.® 
by their bare counterparts, and First note 

the spin-dependent elastic contribution to the pseudo- 
fermion self energy, shown in Fig. [l](c) S/((t) 
a pJzTn sinh^^ {— D / jm) , where p is the density of states 
(DOS), D is the bandwidth, and we assumed a constant 
DOS, A^O: in the band. Since we assume the Fermi surface 
is not nested, opening of the gap at parts of the FS close 
to the magnetic Brillouin Zone boundary leads to a build 
up of density of states, D±{e) for the and '— ' bands re- 
spectively, at the gap edge E±{'k) ^ 7 |m|. The detailed 
shape of the DOS depends on the underlying band struc- 
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FIG. 2: Density of states for the two bands. For all numerical 
calculations, except when noted otherwise, the bandwidth is 
taken to be 10 and 7m = 0.05. 
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FIG. 3: Flow of v{--) and Jz{ ) where the initial values 

are Jz = —0.4 and Jx = —0.4. For comparison we also show 
the flow of Jzp(— )in the absence of AFM order. 



ture of the metal. Here, for simplicity and without loss 
of generality, we model D±{e) ~ Di{e) + D2±{e), where 
Di(e) = iVomin(|e|/7m, 1) is the linearly suppressed DOS 
below the AF gap edge, and the Lorentzian contribution, 
D2±(e) = (a7TO/7r) [(e ^ 7771)^ + (7771)^] . For weak 
AFM, 770. ^ D, the choice a — '-fmNo/2 conserves the 
number of states. The resulting DOS for the particle-hole 
symmetric case, Z?_(— e) = £'-|_(e), is shown in Fig. [51 
most of the lower(upper) band is occupied(empty). De- 
spite this complex behavior of the density of states, eval- 
uating S/(cr) with constant DOS is qualitatively accurate 
since there is no (logarithmic) frequency dependence to 
the self energy at this order. We use the full model den- 
sity of states in computing the renormalization group flow 
and carrying out the variational ansatz below. 
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To one loop we find that T"^ = —^3, while to O(J^) 
0. Both in the particle-hole and particle-particle 
channels the terms cr^S'^ and cr^S'° appear with oppo- 
site signs simply due to spin commutation rules. How- 
ever, if the impurity is at a site with non- vanishing local 
magnetization, an additional relative sign change for the 
£7^5° term between the two diagrams in Fig. [T] is due 
to the splitting of the band fermion spin states at the 
location of the impurity; hence the contributions of the 
two diagrams add in this channel, while canceling in the 
a^S^ channel. To this order in perturbation theory, the 
spin dependent self energy of the pseudofermions is zero, 



and hence F^^ 



-T^; this is not true at higher orders. 



The particle-particle and particle-hole contributions to 
r4 cancel to all orders due to particle-hole symmetry in 
this channel (recall that the magnetization m is chosen 
along z). In the absence of particle-hole symmetry, how- 
ever, this spin non-conserving term is a relevant pertur- 
bation about the ordinary Kondo effect. With particle- 
hole symmetry, F^ and are the only such relevant 
perturbations, while the asymmetry introduced through 
F|" is only marginal at the fixed point. 



D. Renormalization Group (RG) Flow. 

We now consider the renormalization of the coupling 
constants by integrating out the fermion states near the 
band edge. As explained above, the spin non-conserving 
(T^S'° interaction channel appears in the vertex. To in- 
vestigate its flow, we introduce an effective exchange field 
acting only on the fermions: even though such an inter- 
action is absent from the starting Hamiltonian, it is gen- 
erated at one loop level and determines the RG flow by 
polarizing the conduction band at the impurity location. 
To this end we add to the Hamiltonian the term 



27V 



S'°cra^/„kQ/*,k'^a|ji^aa„,i^'^ . (11) 



k,k'. 



In solving the flow equation we impose the initial con- 
dition that V = and follow its evolution. Defining 
h (w) = ^mj \JliJ^ -\- 7^777^, we find the RG equations for 
the coupling constants to order 0{,P ^ w^, Jv) 



2N V (2A^)^ 



/ N vJz , , ,\ Suj , 



(27V)2 



UJ 



SJj_ .JzJ± Suj {v)J± h{uj)6Lj 

(5w 



(27V) 
where, 



.9(w) 



(2iV)2 UJ 
( 



2(2iV)2 
y \ h (lo) 5lo 



to 



\{2NY {2Nf 



(14) 



f{u) = {D+{Lu)+D^{-Lj) + D+i-cu) + D^iuj)), 
g{u;) = {D+{uj) + D^{-u:)-D+{-oj)-D^(u:)), 

and y — -2' ~ Note that for a paramagnetic metal 
with particle-hole symmetry g{uj) = 0, and the equations 
reduce to those for the usual Kondo effeclii. 

The flow of the exchange field v and the coupling con- 
stant Jz is shown in Fig. [3] for the density of states from 
Fig. [21 For comparison we also show the flow of Jz^ in 
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the absence of antiferromagnctic order. Note that both 
the Kondo coupling J, and the vertex v grow, and, at the 
one-loop level (with no self energy corrections), flow to 
strong coupling. This indicates that the Kondo screening 
is modified in the AFM state. To understand the nature 
of the fixed point we perform a calculation using a vari- 
ational wavcfunction motivated by the one loop results. 



IV. VARIATIONAL ANSATZ. 

To estimate the moment, binding energy and the impu- 
rity density of states in the ground state, we diagonalize 
the Hamiltonian within a restricted subspace spanned by 
the states (4 and 11 are the states of the impurity spin, 
and n = ± as before) 



Ink 



2nk 



1 

1 

71 



(wk - nuk) alkT® -U- ~ (^'k + nuk) a\^y® 
(wk + nuk) alkT® -U- + ("k - "-Wk) alki® "fr 



The states chosen here generalize those used to de- 
scribe the ground state in the standard Kondo problen*i^, 
which yield the same result for the ground state as the 
large N approximationi^ or slave Boson methodic. In 
particular, the state ipi for Uk = l,i'k = reduces to 
the singlet formed between the conduction electron and 
the impurity spin, while -02 in the same limit becomes 
the triplet component. In contrast to the paramagnetic 
metal, here these two states are coupled by the exchange 
interaction, and hence both need to be considered when 
constructing the variational wave function. 

We approximate Wk and for all momentum states 
with the same energy, e, by and v^, and work in 
basis states labeled by energy. Diagonalization within 
this subspace gives the eigenstates in the form: -0 = 
X^T-e ''') ^iT-e + ^ ''') ^2Te)- We use the density of 
states of the band fermions. Fig. [21 to numerically evalu- 
ate the energy and the variational coefhcients a and b of 
the eigenstate with the lowest energy. In Fig. [H we plot 
this energy for different values of the parameter 7m and 
bare isotropic exchange coupling, J. 

As expected, the binding energy decreases with in- 
creasing magnetic order m. We also estimate the total 
moment of the collective state by using the fcrmion spin 
operator, s = (1/2) Y.\^,n «ika'^a/3^"k/3, and plotting the 
net moment (S -|- s)^ as a function of m in Fig. 21 Here 
S is the impurity spin operator, and the unscreened case 
corresponds to = 3/4. As the antiferromagnctic order 
increases, the screening is less effective and the moment 
at the impurity site increases. Note significant change 
from the full screening in the absence of AFM order, 
[S + (tY = when m = 0, to a substantial unscreened 
moment at = 0.01. This results from the extreme 
sensitivity of Kondo screening to the low energy features 
of the conduction band. 



0.50- 



^ 0.45 

c 
o> 
E 

O 0.40 

E 
ra 

o 0.35 




— ■— moment 
— ■— energy 



0.00 0.02 0.04 0.06 



FIG. 4: Binding energy normalized to that in the absence of 
magnetic order, Eo{'ym — 0) and net moment of the collective 
state for fixed value of the exchange coupling, J = —0.04. 
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FIG. 5: Local spin-resolved DOS for the impurity located 
at the "up-spin" site for spins parallel(up) and antiparal- 
lel(down) to m where 7m = 0.05 and J = —0.04. 



DENSITY OF STATES. 



To estimate the single particle density of states we 



)^^. The plot of the 



compute n{e,a) = (X]„«^ 

occupation for the two spins is shown in Fig. [5l Since 
the impurity sits at a site with the full symmetry of the 
square lattice, m defines the magnitude and direction of 
the electronic moment at the site in the absence of the 
impurity. Fig. [5l shows that the collective state formed 
by the impurity and the fermions leads to a splitting of 
the impurity density of states. The spin state antiparallel 
to m hybridizes with the fermions to give a rather broad 
peak shifted away from the chemical potential, while for 
the spins parallel to the sublattice magnetization, the 
density of states shows two peaks, one at the chemical 
potential while the other at the antiferromagnctic energy 
scale TO. We expect that the spin-resolved local density 
of states in the antiferromagnctic state shows the same 
salient features even if computed fully sclf-consistcntly. 

Our calculations have been performed in the mag- 
netically ordered phase. We ignored the spin wave 
contribution, which becomes important near the mag- 
netic phase transition; the effect of critical fluctuations 
leads to a multichannel Kondo problem near the crit- 
ical poin1^. Coupling to the spin waves introduces a 
term of order J^gx to one loop where g is the coupling 
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of order jm and x the spin susceptibility of the form 
l/(a;^ + |q — Q|^+C~^)- Since in the ordered phase ^ is fi- 
nite for the modes that couple to the fermions in the long- 
wavelength limit, the contributions due to spin waves are 
higher order in perturbation theory (0(7771 J^^^)), and 
can be safely ignored here. 

VI. SUMMARY. 

To summarize we find that in an itinerant antiferro- 
magnet, the Kondo screening of the impurity moment 
competes with a spin non-conserving coupling, which 
originates from the spin structure of the quasiparticles 
eigenstates in the AFM phase, and the nature of which 
depends on the local symmetry of the impurity site. Our 
result implies that in the heavy fermion AFM state the 
Kondo screening is incomplete; this is in contrast to a 
very recent variational Monte Carlo calculation for the 



Kondo lattice that found screening of the moment across 
the transition from paramegnet to antiferromagnet^. 
Wc have calculated by a simple variational method the 
residual ground state moment and its distribution as a 
function of energy which mimics what occurs in the an- 
tiferromagnetic heavy fermion state. This is the first 
step in understanding the ordered phase of Lattice Kondo 
model and its approach to criticality. 
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